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, Abstract 



The relationship between deexcitation energies of superdeformed secondary minima 
relative to ground states and the density dependence of the symmetry energy is investi- 
' gated for heavy nuclei using the relativistic mean field (RMF) model. It is shown that 

the deexcitation energies of superdeformed secondary minima are sensitive to differences 



in the symmetry energy that are mimicked by the isoscalar-isovector coupling included in 
. the model. With deliberate investigations on a few Hg isotopes that have data of deexci- 

tation energies, we find that the description for the deexcitation energies can be improved 
due to the softening of the symmetry energy. Further, we have investigated deexcitation 
energies of odd-odd heavy nuclei that are nearly independent of pairing correlations, and 
have discussed the possible extraction of the constraint on the density dependence of the 
symmetry energy with the measurement of deexcitation energies of these nuclei. 

PACS: 21.65.Ef, 21.10.-k, 21.60.Jz 
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1 Introduction 

The nuclear symmetry energy plays an important role in astrophysics [1, 2, 3, 4], the structure 
of neutron- or proton-rich nuclei, and the reaction dynamics of heavy-ion collisions, see, e.g., 
Refs. [5, 6, 7]. However, the density dependence of the symmetry energy is still poorly known, 
for instance, see Ref. [7]. Recently, considerable progress has been made in constraining the 
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density dependence of the symmetry energy using data from heavy-ion reactions [8, 9, 10, 11, 
12]. On the other hand, it is promising to constrain the density dependence of the symmetry 
energy at subsaturation densities by accurately measuring the neutron skin thickness in heavy 
nuclei. In the past, Horowitz ct. al. proposed to measure the neutron radius of 208p|-, 
by virtue of the parity-violating electron scattering on the neutrons that promises a 1% 
accuracy [2, 5, 13, 14]. While the precision measurement is still in progress, it is valuable to 
explore whether some other structural properties of finite nuclei are sensitive to differences 
in the symmetry energy. 

As one knows, the neutron skin thickness in heavy nuclei depends sensitively on the density 
dependence of the symmetry energy. The sensitive probe to differences in the symmetry 
energy may thus possibly exist in systems that undergo a relative variation of the proton and 
neutron matter distributions. This relative variation can follow from the collective excitation 
or deexcitation between the ground state and superdeformed secondary minimum (SSM) 
in heavy nuclei with appreciable neutron excesses. Especially, the large relative variation 
can be expected for nuclei in the mass region A ~ 190 where the prolate superdeformation 
of the secondary minima usually occurs with the oblate ground states [15, 16, 17, 18]. In 
deed, quite different deexcitation energies of the SSM for nuclei in the region A ~ 190 were 
predicted by various models that usually diversify the symmetry energies, see Ref. [19] and 
references therein. Since in obtaining the deexcitation energy the isoscalar ingredients of the 
ground state and SSM cancel largely, the variation of the deexcitation energy can be mainly 
attributed to the uncertainty of the symmetry energy that is controlled by the isovector 
potential. However, a direct relationship between the deexcitation energy and the density 
dependence of the symmetry energy is not available in the literature. In this work, it is thus 
meaningful to establish this relationship. It is the aim of this work to constrain the density 
dependence of the symmetry energy in the relativistic mean-field (RMF) model through 
deexcitation energies measured and to be measured for nuclei in the region A ~ 190. We 
note that the constraint on differences in the symmetry energy has been investigated using 
properties of the isovector giant and pigmy dipolc resonances [20, 21, 22, 23]. Apart from 
these dynamical resonances, the SSM features a static structure. Indeed, this is an attempt 
to constrain the density dependence of the symmetry energy using the information of atomic 
masses, since the deexcitation energy is the difference between masses of the SSM and ground 
state. 

The paper is organized as follows. In Section 2, we briefly introduce the formalism of the 
deformed RMF model for finite nuclei. Results on the SSM and ground states of finite nuclei, 
especially the deexcitation energies are presented in Section 3. A summary is finally given in 
Section 4. 
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2 Formalism 



The model lagrangian is written as: 

L = i^ii^^d" - Mat + g^a - ^o^Tm^" - 5p7/.T36^ ~ 4 ^ ^3)7/.^''] 

^\{d^c7d^a - mW) + Uia,u;'', b^), (1) 

where ip, cr,uj, and 60 are the fields of the nucleon, scalar, vector, and neutral isovector- vector, 
with their masses Mjv, rwo-, m,^, and nip, respectively. is the photon field. gi{i = a,uj,p) 
are the corresponding meson-nucleon couplings. Fni,, B^u and Afj^i, are the strength tensors 
of oj and p mesons, and photon, respectively 

F^u = dp,ujy - dyUfj,, B/^^ = - dybon, A^^, = d^Ay - dyA^,. (2) 

The self-interacting terms of cr, a; mesons and the isoscalar-isovector coupling are given gen- 
erally as 

+4glglAyCo^oj''boi,b^. (3) 

Here, the isoscalar-isovector coupling term is introduced to modify the density dependence 
of the symmetry energy. In fact, the symmetry energy can be modified by many theoretical 
factors, for instance, the isoscalar-isovector coupling terms[2], the isovector-scalar mesons, 
the density dependent coupling constants [24], and the model chirality constraint [25]. In the 
past, the density dependence of the symmetry energy had been extensively explored through 
the inclusion of the isoscalar-isovector coupling terms [2, 3, 5, 13] in the RMF theory, and 
this allowed one modify the neutron skin of heavy nuclei without compromising the success 
in reproducing a variety of ground-state properties [2, 5]. In this work, the isoscalar-isovector 
coupling term is thus included in the deformed RMF model to modify the density dependence 
of the symmetry energy. 

Using the Euler-Lagrangian equation, the Dirac equation of motion in RMF is written as 

[-ia ■ V + PM^ + g^uoir) + ^pTa^oCr) + e^{l + r3)^o(r)]V'i(r) = F«Vi(r), (4) 

with = Mjv — ga(^{^) and being the nucleon eigen energy. For simplicity, the isospin 
subscript for the p-meson field is omitted hereafter. For the mesons and photon, the equations 
of motion are given as 

{A-ml)ct^{y) = -s^{v) (5) 
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where for the photon, = 0, and 



ffpP3(r) - Sglgl Ayboir)uj^{r), rho , 



(6) 



epc(r), 



photon. 



Here ps, pb, Ps and pc are the scalar, vector, isovector and charge densities, respectively. The 
total binding energy is written as 



total 



En + Ecr + £'a;o + Ejjg + Ec + Eqm 

Y,{E^ -MN)-\j d\[g,a{v)p,{v) + ^-g2a\v) + \g3a\v)] 



(fr[g^u^Q{v)pB{v) + ya;o(r)] 



gp / ci^r6o(r)[y93(r) + 8c;pC;^AvW^(r)6o(r) 



-\e J dVAo(r>e(r)-^41AV3. 



(7) 



In practical calculations, we also include the BCS pairing interaction using the con- 
stant pairing gaps which are obtained from the prescription of Moller and Nix [26]: A„ = 
4.8/7V^/3, Ap = with N and Z the neutron and proton numbers, respectively. Here, 

the cut-off 82 A~^/^ MeV above the nucleon chemical potentials is used to normalize the pair- 
ing energy. We note that the BCS description for nucleon pairing is not as popular as the 
Bogoliubov approach, e.g., see Refs. [27, 28] and references therein. However, in this work we 
do not pursue a complete description of nucleon pairings but investigate the relative change 
in binding energies of the ground state and SSM with respect to the density dependence 
of the symmetry energy. For both ground states and SSM, the numbers of oscillator shells 
Np = Nb = 18 are used in the basis expansion. The solution of the coupled Dirac and 
meson equations can be obtained in an iterative procedure that can be easily found in the 
literature [29], and it is not reiterated here. 



3 Results and discussions 

In this work, we make analyses based on the calculations with the RMF parameter set 
NL3* [30] that is an improved version of the parameter set NL3 [31]. Also, the isoscalar- 
isovector coupling is included to modify the density dependence of the symmetry energy. This 
new parameter set has been successfully tested by the properties of the ground states and 
collective excitations of fintie nuclei. Based on this parameter set, here we investigate the 
deexcitation energies relative to the ground states for Hg and Au isotopes, while attention will 
be paid to the sensitivity of the deexcitation energy to differences in the symmetry energy. 
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Prior to practical calculations, it is useful to write here the explicit expression of the 
symmetry energy in the RMF models 

where m* is the p-meson effective mass with m* = ^Jfrfip + ^^v{guj9 p^oT ^ is the isospin 
asymmetry with 5 = pz/ Pb, and Ep is the Fermi energy. The first term is the potential part 
of the symmetry energy, and the second term is the kinetic part. The modification to the 
symmetry energy is dictated by the potential part through the isoscalar-isovector coupling. 

Since the symmetry energy is not well constrained at saturation density, some average of 
the symmetry energy at saturation density and the surface energy is needed according to the 
constraint from the binding energy of nuclei. For a given Ay, wc follow Ref. [2] to readjust the 
pNN coupling constant Qp so as to keep the symmetry energy unchanged at kp = 1.15 fm~-^. 
In doing so, the symmetry energy is softened by the isoscalar-isovector coupling, as shown in 
Fig. 1. As a result, it was found that the neutron skin thickness relies sensitively on the Ay, 
while the total binding energy of ^'^^Pb just changes by a few MeV with the Ay of interest [2]. 
The small change in the total binding energy results from two cancellation mechanisms. 
The first mechanism is inherent in the RMF theory due to the cancellation between the big 
scalar attraction provided by the a meson and the big vector repulsion provided by the oj 
meson [32]. Though the isoscalar-isovector coupling modifies predominantly the isovector 
potential, its infiuence can be passed on to the isoscalar potential through modified nucleon 
matter distributions. Here, the a field energy changes with respect to Ay coherently with 
the uj field energy but with an opposite sign. The second mechanism can be understood 
according to the virial theorem, which can be written in the non-relativistic form as: 

where V{r) is the nuclear potential. As shown in Fig. 1, the symmetry energy is modified 
oppositely below and above the fixed point {kp = 1.15fm~^), and this is ensured by the similar 
modification to the isovector potential. Due to the strong coupling between protons and 
neutrons, modifications in the isoscalar potential are similar, albeit small. According to the 
virial theorem, the modification to the nucleon kinetic energy is associated with the gradient 
of the potential. The opposite modifications to the potential give rise to opposite increments 
for the gradient of the potential. Thus, the cancellation between opposite increments leads 
to just a small change to the nucleon kinetic energy. Eventually, the change in the total 
binding energy is small with respect to Ay. Nevertheless, since the empirical binding energies 
of finite nuclei were firstly reproduced quite accurately by the best-fit models, the variation 
of the total binding energy, albeit small, can not simply be used to constrain the density 
dependence of the symmetry energy. 

If the variation of the binding energy is expected to be informative to the density depen- 
dence of the symmetry energy, it should be the relative variation of the binding energies of 
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Figure 1: Symmetry energy as a function of the density with the NL3*. The various Ay are 
denoted for corresponding curves. 

two systems that clearly differ by an isovector ingredient. In this work, we are thus inter- 
ested in the deexcitation energy of the SSM relative to the ground state. However, since the 
isoscalar-isovector couphng modifies not only the isovector potential but also the isoscalar 
potential, an efficient probe to the density dependence of the symmetry energy should be 
able to separate clearly the modification to the isovector potential from that to the isoscalar 
potential. In the following, it is necessary to exhibit such a separation in obtaining the de- 
excitation energy. In order to analyze the potentials obtained with the deformed RMF code 
conveniently in one dimension, it is necessary to perform the multipole expansion as 

U{r) = V4^ J2 UL{r)YLo{e), (10) 

L=0,2,4,- 

where Ylq is the spherical harmonic function. The deexcitation energy, which is the energy 
difference between the secondary minimum (s.m.) and the ground state (g-S.), is associated 
with the corresponding difference of the potentials 

AUr{r) = Ui%^{r)-Ui^^"g,{r), (11) 

where the superscript iso represents the isoscalar (IS) or isovector (IV). Here, the isoscalar 
and isovector potentials are given as Ul^{r) = gujOj{r) — gaO'i'f) and ?7£^(r) = gpbo{r), 
respectively. 

Fig. 2 displays the difference of potentials between the SSM and ground state for I92,i94jjg 
as a function of radius. In Fig. 2, only the two most important components Uo{r) and U2{r) 
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Figure 2: Difference of the potentials between the SSM and ground state of I92,i94jjg 
function of the radius for various Av The calculation is performed within the NL3*. The 
left panel is for the isoscalar potential and the right panel for the isovector one. The result 
for ^^^Hg is displaced downwards by 2 MeV for clarity. 

that are respectively responsible for the spherical matter distribution and quadrupole defor- 
mation are drawn, while the difference between the higher-L components of the potentials is 
becoming clearly smaller and much less sensitive to the Av. It is shown in Fig. 2 that the 
difference of the isoscalar potential is almost independent of the isoscalar-isovector coupling 
Ay. On the other hand, the difference of the isovector potential varies significantly with the 
Av in the surface region. Especially, a large variation against the Av is seen for the isovector 
potential difference AC/2(r), associated with the large oblate-prolate shape difference between 
the ground state and SSM in ^^^Hg and ^^^Hg. These naturally exhibit a clear separation 
of the relative variation of the isovector potentials from that of the isoscalar potentials. In 
particular, the relative variation of isoscalar portions of the two states with respect to the Av 
nearly cancels out. 

Since the isoscalar-isovector coupling modifies just slightly the nucleon potentials, it can 
be treated as a small residue interaction. This rules out the occurrence of large coherent 
changes in nuclear potentials of the ground and metastable states while with slight relative 
changes between various states, and thus the small change in total binding energies can be 
determined by that of nuclear potentials with rather stable nuclear structures. Consequently, 
the shift of the deexcitation energy caused by the isoscalar-isovector coupling can be obtained 
from the small relative change in nuclear potentials between the ground state and SSM. Due 
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to the nearly exact cancellation between the variations of isoscalar potentials in the SSM 

and ground state, the variation of the deexcitation energy results predominantly from the 
modification to the isovector potential caused by the isoscalar-isovector coupling. Besides 
from the potential part, the kinetic energy and the nonlinear term of the isoscalar-isovector 
coupling also directly cause the variation of the deexcitation energy, but they can after all be 
determined from the nuclear potential, for instance, according to Eqs.(5) and (9). Eventually, 
this builds up a direct relationship between the uncertainty of the deexcitation energy and 
the modification to the symmetry energy. 

Table 1: Properties of the ground state and SSM for ^^^Hg within the NL3* with respect 
to Av. The binding energy per nucleon (B/A), charge radius (rc), neutron skin thickness 
{fp — fn)i and quadrupole deformation parameter j3 are listed. The properties of the SSM 
are denoted by the subscript SD. Energies are in unit of MeV and radii in unit of fm. 



Av 


9p 


B/A 


rc 






(B/A)5D 


{rn - rp)sD 


f^SD 


0.000 


4.5748 


7.913 


5.460 


0.213 


-0.146 


7.885 


0.182 


0.625 


0.010 


4.9005 


7.928 


5.461 


0.189 


-0.145 


7.898 


0.159 


0.622 


0.020 


5.3074 


7.940 


5.464 


0.165 


-0.144 


7.907 


0.137 


0.618 


0.030 


5.8360 


7.948 


5.468 


0.141 


-0.142 


7.914 


0.115 


0.614 



In Table 1, we tabulate as an example calculated quantities of the ground state and 
SSM of ^^^Hg with respect to the isoscalar-isovector coupling constant Ay. We see that 
except for the neutron skin thickness all the properties of ^^^Hg depend just slightly on 
the Ay. With the inclusion of the isoscalar-isovector coupling, the neutron skin thickness 
reduces appreciably. These arc consistent with the early findings for heavy nuclei in the 
literature [2, 5]. It is seen in Table 1 that the SSM has a different neutron skin thickness from 
that of the ground state. This is due to the excitation of protons and neutrons in the vicinity 
of Fermi surfaces to different shells. Specifically, in ^^^Hg neutrons with the configuration 
2/|y2 1^9/2 t'^^t below the major shell N = 126 in the ground state are excited to intruder 
orbitals that originate from spherical orbitals 2^9/2 liii/2 above the major shell A'" = 126. The 
neutron occupation in intruder orbitals forms the important configuration '^gg/2^'''ii/2 
superdeformation in the SSM. For protons, similar excitation to intruder orbitals occurs from 
the shell 50 < Z < 82 in the ground state up to the shell Z > 82 in the SSM. Besides the 
characteristic excitation between the major shells, the excitation in subshclls also contributes 
to the formation of the superdeformation in the SSAI. Thus, the creation of interior holes 
and occupation of exterior orbitals at different major shells for protons and neutrons results 
in an appreciably different neutron skin thickness in the SSM from that in the ground state. 
As a result, the SSM differs from the ground state by an isovector ingredient from which the 
sensitivity of the deexcitation energy to the isoscalar-isovector coupling originates. For other 
Hg isotopes, the results are similar to those for ^^^Hg, and are not given here. 

In Table 2, we list the deexcitation energies relative to the respective ground state for 
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Table 2: Deexcitation energies (MeV) of the SSM relative to ground states for a few Hg 
isotopes with respect to the Ay. The experimental values for I92,i94jjg ^^j^^^ isiRg are taken 
from Refs. [19, 33], respectively. 



Model 


Av 




^^%g 


I'^lHg 


NL3* 


0.000 


5.52 


4.03 


3.59 




0.010 


5.81 


4.40 


3.94 




0.020 


6.27 


4.73 


4.26 




0.030 


6.60 


5.03 


AAA 


Expt. 




6.0 


5.3 


4.6 



I94jjg^ 192 jjg g^^^ 191 jjg ii^Q^^ have experimental data [19, 33]. We may observe from Table 2 
and Fig. 1 that the deexcitation energies are sensitive to differences in the symmetry energy. 
Comparing with the experimental data, it is seen that the description for the deexcitation 
energies can be largely improved by including the isoscalar-isovector coupling. With Av = 
0.03, the deexcitation energies for ^^^Hg, ^^^Hg and ^^^Hg agree fairly well with the data 
within acceptable accuracy. This confirms the softening of the symmetry energy. 

Table 3: Deexcitation energies for Hg isotopes recalculated with the slightly modified rua 
with respect to the Ay. As compared to gp in Table 1, it is just slightly modified to keep the 
symmetry energy unchanged at kp = 1.15fm^^. Ground-state binding energies per nucleon 
B/A are hsted for ^^^Hg. 



Av 


rria (MeV) 


9p 


194Hg 


192jjg 


191jjg 


B/A 


0.000 


502.5742 


4.5748 


5.51 


4.06 


3.59 


7.913 


0.010 


502.6050 


4.9006 


5.63 


4.37 


3.85 


7.912 


0.020 


502.6310 


5.3076 


6.19 


4.62 


4.29 


7.913 


0.030 


502.6465 


5.8363 


6.54 


4.99 


4.46 


7.913 



Now, we clarify the concern whether the variation of deexcitation energies can be dis- 
guised by the much larger change in the ground-state energy, which is about 7 MeV for the 
Av of interest, as estimated in Table 1. To do this, we just need to observe the variation 
of the deexcitation energy in the case of the constant ground-state energy. It is possible to 
reduce or even eliminate the variation of the ground-state binding energy by suitably choos- 
ing the fixed point which can affect the extent of cancellations, whereas this is numerically 
complicated. Alternatively, the elimination of the variation in the ground-state energy can 
actually be fulfilled by slightly readjusting the meson-nucleon coupling constants or meson 
masses. Without priority, here we realize it by slightly readjusting the a meson mass rria- 
The recalculated results are given in Table 3. As seen in Table 3, the variation of the total 
ground-state binding energy is nearly eliminated by just modifying the mg- up to 0.07 MeV. 
With this slight readjustment of the rua, the modification to the incompressibility is just 
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about 0.8 MeV. Noticeably, the variation of deexcitation energies is almost unchanged, as 
compared with that given in Table 2. Definitely, the sensitivity of deexcitation energies to 
differences in the symmetry energy is irrelevant to the variation of the ground-state energy. 
The underlying physics for this is attributed to the fact that the variation of the deexcita- 
tion energy is conditioned predominantly on the modification to the isovector potential that 
changes the density dependence of the symmetry energy, as analyzed for results shown in 
Fig. 2. 

On the other hand, the difference between deexcitation energies of ^^^Hg and ^^^Hg is still 
large (about 1.5 MeV), as compared to the experimental value 0.7 MeV, and is not reduced by 
including the isoscalar-isovcctor coupling. It means that the accurate deexcitation energies 
for ^^^Hg and ^^^Hg both can not be obtained with the same It is seen in Table 2 
that the difference between the deexcitation energies for ^^^Hg and ^^^Hg is not sensitive 
to the Ay. This is because the isospin asymmetries 5 in ^^^Hg and ^^^Hg are not very 
different. Consequently, the considerable reduction of the difference between the deexcitation 
energies for ^^^Hg and ^^^Hg should resort to the alteration of the isoscalar potential. For 
instance, this difference is reduced to be 0.5 MeV with the RMF parameter set TMl [34] whose 
isoscalar potential is moderately different from that with the NL3* due to the inclusion of 
the a;-meson self-interaction. However, with the w-meson self-interaction, the softened vector 
potential leads to small deexcitation energies and much shallower potential wells of the SSM. 
As deexcitation energies for ^^^Hg and ^^^Hg are both much smaller than the data within the 
TMl, the isoscalar-isovector coupling can be included to partly compensate this discrepancy, 
whereas an overall compensation seems yet to be available in the model that features the 
cj-meson self-interaction. 

Further, we notice that the deexcitation energy is nearly of the same magnitude as the 
pairing energy. The correct deexcitation energy of some nuclei may be obtained with the 
inclusion of the isospin-dependent pairing interactions. For instance, using different pairing 
gaps and strengths in Ref. [35], the deexcitation energy of ^^^Hg with the NL3 was repro- 
duced in nice agreement with the experimental value. On the other hand, the difference of 
the deexcitation energies between ^^^Hg and ^^^Hg was still 0.8 MeV larger than the experi- 
mental value [35]. If this prescription is applied to resolve the discrepancy of the deexcitation 
energy with the experiments for other isotopes, it will impose stringent constraints on pair- 
ing strengths. On the other hand, the symmetry energy in finite nuclei can successfully be 
extracted regardless of the pairing interaction [36] . This implies that the observable response 
to the symmetry energy would be rather independent of the pairing interaction. Specifically, 
we find that the deexcitation energies of the SSM in odd-odd nuclei in the region A ~ 190 
just depend very weakly on the pairing interaction, while the sensitive dependence on the 
isoscalar-isovector coupling does not alter at all. The former is because the pairing interac- 
tion in odd-odd nuclei is substantially suppressed by the Pauli blocking and the nucleonic 
current that breaks the double degeneracy. If the pairing interaction turns off while reviving 
the nucleonic current and the contribution of the tt meson in these odd-odd nuclei [37], the 
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Figure 3: Deexcitation energies (i?-^) of the SSM for odd-odd Au isotopes versus the neutron 
skin thickness in ^°^Pb with the NL3*. The results are obtained by varying the Ay properly 
with Qp. The mass number is marked for each curve. 

deexcitation energy is just changed by about 0.2 MeV. Therefore, no appreciable effect on 
deexcitation energies can be observed by modifying pairing strengths in these odd-odd nuclei. 
In Fig. 3, wc display as examples the deexcitation energies for a few odd-odd Au isotopes as a 
function of the neutron skin thickness for ^*^^Pb. Since the neutron skin thickness in ^OSp]-, ^jg, 
pends almost linearly on the Ay [2, 5], the sensitive dependence on the Ay can be observed for 
deexcitation energies of odd-odd Au isotopes. We can thus expect that the measurement of 
the deexcitation energy of these odd-odd isotopes can be useful not only for constraining the 
density dependence of the symmetry energy but also clarifying this theoretical uncertainty. 

Next, let us discuss more the factors that are necessary for the description of the deexci- 
tation energy. First of all, since the SSM is a result of the collective excitation that undergoes 
a significant isovcctor change, the appropriate isovector potential of the model is necessary. 
We find from Table 2 and 3 that the description for deexcitation energies of various isotopes 
is improved by including the isoscalar-isovector coupling. In other words, the inclusion of 
the isoscalar-isovector coupling improves the isovector potential of the model. Secondly, we 
note that the isoscalar potential is also an important ingredient to carry out correct deexcita- 
tion energies. Especially, the difference between the deexcitation energies for the neighboring 
even-even isotopes is tightly associated with the appropriate fit to the isoscalar potential. The 
present study indicates that an overall description for the deexcitation energies and difference 
between them needs to reconstruct the model appropriately with the isovector potential that 
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brings out a softened symmetry energy. 

At last, we should stress that though the present work is performed based on a specific 
parametrization (NL3*), the conclusion that the deexcitation energy is sensitive to differences 
in the symmetry energy is rather free of this specific RMF model. We have examined this 
with many other RMF models. The current situation is that most best-fit models have been 
constructed regardless of the detail of the density dependence of the symmetry energy. Thus, 
the agreement with data of deexcitation energies may play a significant role in reconstructing 
models with the appropriate density dependence of the symmetry energy. Moreover, we 
have noted that the treatment of the nucleon pairing with the BCS theory is comparatively 
simple as compared with the Bogoliubov approach. This would affect the accuracy of the 
correlation between deexcitation energies of even miclci and the density dependence of the 
symmetry energy. One can expect the improvement with the relativistic Hartrce-Bogoliubov 
models. However, this is beyond the capacity and scope of this preliminary effort. 

4 Summary 

In summary, we have studied the relationship between the deexcitation energies of the SSM 

in heavy nuclei and the density dependence of the symmetry energy. The investigation is 
based on the RMF model (NL3*) with the isoscalar-isovcctor coupling included to modify 
the density dependence of the symmetry energy. It is found that the uncertainty of the de- 
excitation energy originates almost uniquely from the modification to the isovector potential 
induced by the isoscalar-isovcctor coupling. The deexcitation energies can thus serve as a 
theoretical probe to the density dependence of the symmetry energy. As a result, we find 
that the theoretical estimates of the deexcitation energies of Hg isotopes can be improved 
by the inclusion of the isoscalar-isovcctor coupling that softens the symmetry energy. To 
separate the effect of the symmetry energy on deexcitation energies from that of pairing cor- 
relations, we propose to measure the deexcitation energies of odd-odd heavy nuclei in the 
region A ~ 190 such as Au isotopes that are nearly independent of pairing correlations. This 
can be used to constrain the isovector content of the model and thus the density dependence 
of the symmetry energy. 
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